Implications of the Super-K atmospheric data for the mixing angles 6 Vi and #23 
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A three-neutrino analysis of oscillation data is performed using the recent, more finely binned 
Super-K oscillation data, together with the CHOOZ, K2K, and MINOS data. The solar parameters, 
A21 and #12, are fixed from a recent analysis, and A32, #13, and #23 are varied. We utilize the full 
three-neutrino oscillation probability and an exact treatment of the Earth's MSW effect with a 
castle- wall density. By including terms linear in #13 and e := #23 — 7r/4, we find asymmetric errors 
for these parameters #13 = — 0.07ig and e = 0.03^^. For #13, we see that the lower bound is 
primarily set by the CHOOZ experiment while the upper bound is determined by the low energy 
e-like events in the Super-K atmospheric data. We find that the parameters #13 and e are correlated 
— the preferred negative value of #13 permits the preferred value of #23 to be in the second octant, 
and the true value of #13 affects the allowed region for #23. 
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I. INTRODUCTION 



The experimental observation of neutrino oscillations 
implies that at least two of three neutrinos have mass, 
and the mass eigenstates differ from the flavor eigen- 




including the null result from the CHOOZ reactor ex- 
periment [2l[, can be globally understood in terms of 
three mixing angles 0jk, with j = 1,2,3 and j < k, 
one phase 5, and two independent mass-squared differ- 
ences Akj ■= m 2 — to 2 . The separation between two 
of the mass-squared differences is sufficiently large so 
that the data from a given experiment, which may span 
some range of baselines and neutrino energies, can be 
approximately understood within the context of an ef- 
fective two-flavor theory. Experiments detecting solar 
neutrinos [H, 0, i, i, i, @, 0, 0, ©, El and the longbase- 
line (LBL) reactor experiment KamLAND [n], [l2j are 
particularly sensitive to the mixing angle O12 and the 
mass-squared difference A21 assuming the standard rep- 
resentation of the neutrino mixing matrix [22T |. A three 
neutrino analysis [23[ gives a value for the mixing an- 
gle sin 2 6*12 = 0.3041q q34 (2 a error), with a precision 
of 8% at 3(7. The solar mass-squared difference is de- 
termined predominantly by the SNO data [Io| and is 
found to be A 21 = 7.65t° Z x 10 ~ 5 ey2 - Atmospheric 
and accelerator beam stop neutrinos provide experimen- 
talists with a good source with which to measure 623 
and A 32 . MINOS [HI, Gil predominantly determines A 32 
while the mixing angle #23 is determined mainly by the 
Super-K atmospheric da ta [l3l [Til HH [lr| . Present val- 
ues for these parameters [23| are A 32 



eV 2 and sin 2 2 3 = 0.50 
gle, 6>i3 
sin 2 6<i3 
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The remaining mixing an- 
mixes the two scales. This same analysis gives 
< 0.040; recent analyses hint at a value of #13 
differing from zero 



cent anaiys' 
[2IJ2IJ21I 



can be found at RefsTgJlig 
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As we enter the era of precision measurements, global 
analyzes of neutrino data must employ a full three- 
neutrino framework in order to correctly assess the neu- 
trino mixing parameters. This will become evident herein 
as we consider various experiments' impact upon the 
small parameters 613 and e := #23 — 7r/4, the deviation 
of 023 from maximal mixing. The quantitative knowl- 
edge of 6 13 is a particularly important part of neutrino 
oscillation phenomenology because it sets the magnitude 
of possible CP violating effects as well as the size of ef- 
fects that might be used to determine the neutrino mass 
hierarchy. There are presently three new reactor experi- 
ments planned or under construction which are designed 
to measure 6q 3 , Daya Bay [11], Double CHOOZ [H, and 
RENO 3l|; an LBL experiment is also planned, T2K 
f32l ]. The subsequent generation of experiments, which 
will be designed to ascertain the level of CP violation, 
cannot proceed until the current generation better de- 
termines the value of #13. In addition, a more quanti- 
tative knowledge of the mixing angles, and particularly 
of 9i3, can help discern between models and symmetries 
of the physics that underlies neutrino mixing 33] . The 
deviation of #23 from maximal mixing is also important 
in model building as it might indicate the presence of 
a broken symmetry. At short baselines, the oscillation 
probabilities which might probe the mixing angle #13 are 
quadratic in this small parameter; however, we have pre- 
viously shown there are terms in the oscillation probabil- 
ity which are linear and appreciable at very long baselines 
(VLBL) [H, |H, HI] and arise from interference between 
the oscillations driven by the two mass-squared differ- 
ences. This is also a region of the parameter space where 
one can look [3?], [38| for CP violating effects. The sub- 
GeV data set of the Super-K atmospheric experiments 
is potentially sensitive to such effects. Furthermore, it 
was shown in Ref. [36j that there is a non-trivial rela- 
tion between e and 613 for sub-GeV neutrinos at VLBLs. 
As such, the extraction of these parameters from the at- 
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mospheric data requires a full three neutrino treatment 
since approximations overly simplify the correlations of 
the parameters. 

We here investigate atmospheric neutrino oscillations 
with the full three neutrino oscillation probabilities. Be- 
cause we do not use truncated expansions, all terms lin- 
ear in #13 and e will be considered as well as higher or- 
der contributions. We do not expect a large change in 
the extracted parameters as only a limited number of the 
Super-K data bins lie in the region where linear terms will 
be significant. On the other hand, in the context of at- 
mospheric data, #13 is itself a small effect as is the octant 
of #23- Small effects can sometimes have a proportionally 
larger impact on something that is inherently small. In 
keeping with the use of the full three neutrino oscilla- 
tion probabilities, we also utilize the method proposed in 
Refs. [H, E| to treat the MSW effect [HI, E|. With a 
castle wall density profile of the earth, this treatment of 
the MSW effect is exact so that approximate expressions 
for the oscillation probabilities are not needed. We also 
include a model for the multi-ring events, a data subset 
often neglected. 



II. ANALYSIS 

In vacuo, the probability that a neutrino of flavor a 
and energy E v will be detected as a neutrino of flavor j3 
after traveling a distance L is given by 

3 

V a p(L/E u ) = 8 a p - 4 ^2 [UajUakUpkUpj) sin 2 ip jk , 
*=<j, 

j,k=l 

(1) 

with (fjk := 1.27 Aj/j L/_E„, where L is measured in kilo- 
meters, E v in GeV, and the mass eigenvalues in eV. 
The matrix U a i is the unitary matrix that relates the 
mass basis i to the flavor basis a. We assume CP con- 
servation so that the U a i are real. Neutrinos which prop- 
agate long distances through matter of sufficient densi- 
ties can incur significant interactions which are diagonal 
in flavor. For matter of constant density, the upshot of 
these interactions is a modification of the effective mixing 
angles and mass-squared differences so that an oscillation 
formula similar to Eq. ^ holds. The density of the earth 
may be approximated as piecewise constant (43J. In ad- 
dition, for certain energies and densities, the neutrinos 
can undergo parametric resonances in regions of varying 
densities [44(. To account for these interactions, we em- 
ploy a simple model of the earth: a mantle of density 4.5 
gm/cm 3 and a core of density 11.5 gm/cm 3 with radius 
3486 km. Using the methods in Refs. [1^, we are 
able to fully incorporate an exact three neutrino model 
of the neutrino-matter interactions which automatically 
incorporates any possible parametric resonances. 

Our interest is to study and extract the following pa- 
rameters from the experimental data: #13, 823, and A32. 
As such, we fix the solar mixing parameters from a re- 



cent analysis (27|, #12 = 0.58 and A 2 i = 8.0 x 10 -5 eV 2 . 
Given that there is no evidence to indicate CP violation, 
we assume CP conservation and take S = 0. We include 
the details of our analyses of the relevant experiments in 
the Appendix. We comment on them briefly here. The 
Super-K atmospheric data is statistically the most signif- 
icant data set, and it covers a range of over four orders of 
magnitude in L/E. Our analysis employs the most recent 
more finely binned data [l6| , a necessity for studying the 
small parameters #13 and e. We also include the most 
recent MINOS results [3, the K2K results [H, and the 
CHOOZ results [H. 

In order to ascertain the importance of the linear 
and higher-order terms in #13 (and also e), we com- 
pare our results with those generated by the often used 
sub-dominant approximation which arises from an ex- 
pansion in the ratio of the mass-squared differences, 
a = A12/A32. In this approximation, the leading order 
oscillation probabilities are given by 

V ee = 1- Sin 2 2013 sin 2 ((^32) 

V efl = sin 2 6*23 sin 2 26»i 3 sin 2 (ip 32 ) 

V m = 1 - 4 cos 2 6>i3 sin 2 9 23 (1 - cos 2 9 13 sin 2 23 ) 

xsin 2 (v?3 2 ) . (2) 

Additional correction terms [13, [H, [45[ can then be 
added. The results for the sub-dominant approximation 
when compared to the results for the full three-neutrino 
oscillation probabilities will inform us of the size of the 
correction terms. 

We begin by examining the mass-squared difference 
A 32 . We plot Ax 2 versus A 32 using both the sub- 
dominant approximation, Fig. la, and the full three- 
neutrino calculation, Fig. lb, with #13 and #23 as var- 
ied parameters. The [black] solid curves are obtained 
from the Super-K atmospheric data alone. The [red] 
dash-dot curves employ the K2K, MINOS, and CHOOZ 
data, omitting the Super-K atmospheric data. These 
curves are largely determined by the recent MINOS data 
which constrain the mass-squared difference more so than 
Super-K, as is well known. The analysis utilizing all of 
the data sets (atmospheric, K2K, MINOS, and CHOOZ) 
is depicted by the [blue] dashed curves. Notice that 
although the Super-K atmospheric data is not as con- 
straining as MINOS, it combines with MINOS to produce 
a reduced bound, particularly from above. In compar- 
ing the approximation, Fig. la, with the full calculation, 
Fig. lb, we see that the sub-dominant approximation is 
useful for determining the mass-squared difference A32 . 
A very careful inspection will reveal that the full three- 
neutrino analysis produces a slightly larger bound than 
does the sub-dominant approximation. Our results are 
A 32 = 0.25±g;g§ eV 2 at the 90% confidence level. (The 
errors quoted for our calculations will be for A% 2 = 6.25, 
the 90% confidence level for a three parameter fit.) 

We next present A% 2 versus e = $23 — tt/4 using the 
sub-dominant approximation, Fig. 2a, and the full three- 
neutrino calculation, Fig. 2b, with #13 and A32 as var- 
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A 32 (10" 3 eV 2 ) 

FIG. 1: [color online] A\ 2 versus mass- squared difference A32 
for the a.) sub-dominant approximation and b.) full three- 
neutrino calculation. The [black] solid curves utilize only at- 
mospheric data; the [red] dot-dash curves utilize K2K, MI- 
NOS, and CHOOZ data; the dashed [blue] curves utilize all 
the data sets: atmospheric, K2K, MINOS, and CHOOZ. 
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FIG. 2: [color online] A\ 2 versus e for the a.) sub-dominant 
approximation and b.) full three-neutrino calculation. The 
[black] solid curve utilizes only atmospheric data; the dashed 
[blue] curve utilizes all the data sets: atmospheric, K2K, MI- 
NOS, and CHOOZ. 




FIG. 3: [color online] The same as Fig. [T] except Ax 2 versus 
#13 is presented. 



ied parameters. We express our result in terms of 623, 
rather than sin 2 2 023 or sin 2 #23, because the oscillation 
probabilities truly are a function of 6*23. The [black] 
solid curve again represents the Super-K atmospheric 
data alone. The [blue] dashed curve represents the re- 
sults from all data sets: Super-K atmospheric, K2K, MI- 
NOS, and CHOOZ. Adding K2K, MINOS, and CHOOZ 
hardly alters the Super-K result. We do not present the 
results for K2K, MINOS, and CHOOZ alone because this 
data does not yield a reasonable constraint on #23 when 
treated as a linear variable with a varied #13 included in 
the analysis. Only in a two neutrino analysis does K2K 
and MINOS restrict the appropriate variable sin 2 26>23- 
Comparing the sub-dominant approximation with the 
full calculation, we see that the full three-neutrino prob- 
abilities produce an allowed region which is much more 
asymmetric about e = 0. In fact, we find a statistically 
insignificant indication that 023 is greater than tt/4, max- 
imal mixing. The ability of atmospheric data to deter- 
mine the octant of 6*23 has also been investigated in . 
Wc find the value at 90% CL is e = 0.03±g;?§. 

In Fig. 3a, we present A% 2 versus #13 calculated in 
the sub-dominant approximation and full three-neutrino 
formulation with #23 and A32 as varied parameters. Pre- 
viously, it has been shown [161 ] that in the sub-dominant 
approximation the atmospheric data alone restrict #13. 
Focusing upon our sub-dominant calculation, Fig. 3a, the 
[black] solid curve depicts the corresponding result from 
our analysis. As noted in the Appendix, our analysis 
quantitatively reproduces the results in Ref. [16j . For 
A% 2 < 4.6, the 90% confidence level for a two-neutrino 
analysis, we both find sin 2 6»i3 < 0.14 (or |0 13 | < 0.38). 
This is a very important calibration of our analysis tool. 
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The effect of 6*13 on atmospheric oscillations is small, 
and obtaining the same result implies we are reproducing 
small effects, not just the global features of the analysis. 
The dash-dot [red] curve in Fig. 3a is the result of ana- 
lyzing the K2K, MINOS, and CHOOZ data, neglecting 
the Super-K atmospheric data. This curve is mainly de- 
termined by the CHOOZ data. We see that CHOOZ 
is more constraining on #13 than is the Super-K atmo- 
spheric data. However, the dashed [blue] curve presents 
the results utilizing all of the data sets and shows that 
the Super-K data does somewhat reduce the error on 
#13; this is due to the indirect effect arising from Super- 
K further constraining the mass-squared difference A23. 
To obtain the constraints on #13 implied by the Super-K 
atmospheric data, it is important to use the data from 
Ref. [161 ] which is more finely binned than earlier Super-K 
[lH work. Note that the curves in the sub-dominant ap- 
proximation are symmetric about #13 = as is manifest 
from the approximate oscillation formulae, Eq. (T5|). In 
Refs. [H, [23, HfJ > it has been observed that recent data 
imply a statistically insignificant non-zero value for (9 13 ; 
our results are likewise consistent. 

Turning to the full three-neutrino calculation, Fig. 3b, 
we find Ax 2 to be very asymmetric with a strong pref- 
erence for negative #13 when using only Super-K data, 
the [black] solid curve. The [red] dash-dot curve employs 
only K2K, MINOS, and CHOOZ data; it is symmetric 
about the origin so that the asymmetry present when all 
data is included, the [blue] dashed curve, is due to the 
Super-K data. What is more, we see the novel result [47} 
that #13 is constrained from above by the Super-K at- 
mospheric data, not by CHOOZ, while it is constrained 
from below primarily by CHOOZ. 

This conclusion is further reinforced by looking at the 
allowed region for the parameters #13 and #23 as depicted 
in Fig.H We plot the 90% confidence level of Ax 2 = 4.61 
for a two parameter analysis as we fix the third parame- 
ter A32 in calculating these curves. The dash-dot [green] 
curve depicts the results for Super-K atmospheric data 
alone in the sub-dominant approximation. We compare 
this with the dashed [red] curve which also utilizes only 
the Super-K atmospheric data alone but incorporates the 
full three-neutrino probabilities. Again, we see the sig- 
nificant change brought about by incorporating the lin- 
ear, and higher-order, terms in #13. The allowed region 
grows, favoring negative #13. The dash-dot-dot [blue] 
curve utilizes all the data in the sub-dominant approxi- 
mation. It is similar to the dash-dot [green] curve because 
the mass-squared difference A32 is fixed when calculating 
the curves; the main effect of the MINOS experiment is 
to restrict A32. Finally, the solid [black] curve utilizes 
all the data and the full three-neutrino oscillation prob- 
abilities. Note that the upper bound on $13 is similar to 
that from the dashed [red] curve, i.e. the curve also uti- 
lizing the full three-neutrino oscillation probabilities but 
only the Super-K data. For the lower bound on #13 , how- 
ever, we find similarities to the dash-dot-dot [blue] curve 
where the restriction on #13 originates primarily from the 
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FIG. 4: [color online] The allowed region at 90% CL for the 
parameters #13 and 623. The solid [black] curve uses the full 
three-neutrino oscillation probabilities and all the data sets: 
CHOOZ, K2K, MINOS, and Super-K atmospheric. The dash- 
dot-dot [blue] curve represents the use of the subdominant 
approximation and all the data sets. The dashed [red] curve 
depicts the use of the full three-neutrino oscillation probabil- 
ities and only the Super-K data. The dash-dot [green] curve 
depicts the use of the subdominant approximation and the 
Super-K data alone. The vertical (horizontal) straight line 
marks #13 = 0.0 (#23 = 7! "/4). The [violet] star indicates the 
location of the minimum for the analysis that includes all 
the data sets and the use of the full three-neutrino oscillation 
probabilities. 



CHOOZ experiment. Thus we again see that the upper 
bound on $13 no longer arises from the CHOOZ experi- 
ment but is determined by the Super-K atmospheric ex- 
periment, while the lower bound continues to come from 
the CHOOZ experiment. Our final result for this mixing 
angle is 13 = -0.07±°;^. 

The principle effect of utilizing the full three neutrino 
oscillation probabilities is the alteration of the shape of 
the allowed region for #13, particularly the introduction 
of the asymmetry about zero. The absolute minimum 
for x 2 is lowered by only 1.3 [ItJ . This is because the 
minima are very close to #13 = where the linear and 
higher order terms contribute little. 



III. DISCUSSION 

The most striking differences between the sub- 
dominant approximation and the full three-neutrino 
probabilities were seen in the determination of the mix- 
ing angle 6*13, Fig. [3] Additionally, the deviation of 623 
from maximal mixing also produced noticeable features, 
though less striking, Fig. [2] Clearly, the two features are 
non-trivially linked as demonstrated in the allowed re- 
gions depicted in Fig. 0] In fact, we see from Fig. [3] that 
the Super-K data is the source of the asymmetry about 
#13 — in the full three-neutrino model. To flesh out 
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FIG. 5: [color online] The two lowest eigenvalues of the effec- 
tive mass matrix in a constant density mantle as a function of 
neutrino energy E„. The [red] solid curve is (mi) 2 (E v ) , and 
the [blue] dashed curve is (m^) 2 (E v ). 



which subset of the Super-K data results in these asym- 
metries, we examine the various contributions of the data 
to x 2 for a fixed positive and negative value of the mixing 
angle 9 13 taken to be ±0.15. The total difference in Ax 2 
for 9 13 = +0.15 and 6> 13 = —0.15 is ~ 7.0. Focusing on 
the fully contained events, we find that two thirds of this 
change in A% 2 between the positive and negative values 
of the mixing angle comes from the sub-GeV electron-like 
events. Half of the total change in A% 2 (3.5) arises from 
a single angular bin within this subset of data, namely 
the bin for e-like events in which the detected charged 
lepton has zenith angle i? satisfying —0.8 < cosi? < —0.6 
and momentum less than 250 MeV. The detected leptons 
in this bin are produced by neutrinos which travel along 
a very long baseline upward through the earth. Such 
neutrinos fall into the region of L/E where we have pre- 
viously shown effects linear in 6*13 to be sig nificant in the 
oscillation probabilities V eil and |34l . [35l . [36[ ; such 
effects occur in a region where the sub-dominant approx- 
imation is the leading term in an expansion which is not 
convergent. Terms linear in #13 can be even more signif- 
icant [481 ] should an atmospheric oscillation experiment 
be able to take data at energies below 100 MeV. We com- 
ment that the angular bin, — 1.0 < cos i9 < —0.8, is also in 
this L/E region for the low-energy neutrinos; however, a 
typical neutrino which produces leptons in this bin passes 
through the earth's higher density core. We find that the 
core suppresses the oscillations and thus this angular bin 
is not as sensitive to the effects linear in #13 . 

The preference of the data for negative 6*13 can be 
linked to the excess of e-like events in the sub-GeV data 



set |37j . This excess is not present in the /i-like data or in 
the multi-GeV data so that an overall renormalization of 
the atmospheric flux cannot account for the excess. To 
understand the role of the data in extracting #13 and 6*23, 
we examine the relevant oscillation probabilities in the 



limit of a constant density mantle and sub-GeV neutrino 
energies, keeping only terms linear in #13 and e and av- 
eraging over the A32 oscillations; these approximations 
have been discussed previously [H, H3] • As detailed in 
the Appendix, the electron-like events at the Super-K de- 
tector are related to the v e survival probability and the 
conversion probability via lZ e = V ee + f"P ell where r 
is the ratio of the to v e flux at the source. This yields 
the approximate expression 



1l e 



1 



r sin 2 W\ 



I-I + cot^™)^ 
2 r 



sax 



V21- I 3 ) 



Here, 9^2 is the effective mixing angle in matter; addition- 
ally, the phase ip^i employs the effective mass-squared 
difference in matter corresponding to A21. In this ap- 
proximation, we can understand how to effect an excess 
of electron-like events for sub-GeV neutrinos over a long 
baseline, 



i_i + cot(20S; 



13 



e > 0. 



(4) 



Using the same approximations, we can simply express 
the MSW resonant energy 



Er 



A 2 i cos 26*12 
2Vcos 2 6 ) 



13 



(5) 



with V « 1.7 x 10~ 13 eV in the mantle; this yields En on 
the order of 100 MeV. This resonance is apparent when 
we plot the eigenvalues of the effective mass-squared ma- 
trix in the mantle, Fig. the "resonance" is indicated 
by the slight bowing in the curves toward each other and 
is located at the point where the effective mass-squared 
difference is minimal. 

At the resonant energy, one has 6>^ = 7r/4; for neutrino 
energies above the resonance, the effective mixing angle 
in matter increases up to 7r/2. As a consequence, for neu- 
trino energies above 100 MeV in the mantle, the function 
cot(26']5) is negative; to reiterate, the coefficient of the 
6*i3 term in the inequality, Eq. ([4|, is negative. We note 
that for these low energy atmospheric neutrinos r ~ 2 so 
that the first two terms of the inequality approximately 
sum to zero. If 6*13 is restricted to positive values, then 
the mixing angle 6*23 must lie in the first octant (e < 0) 
in order to account for the excess in lZ e . However, if we 
allow 9i3 to run the full range of allowed parameter space 
in a CP conserving theory, then a negative value of this 
mixing angle can easily accommodate the excess in lZ e , 
even permitting $23 to lie in the second octant as is the 
case in our analysis. 

To demonstrate the point regarding the effect of terms 
linear in 6*13, we plot lZ e in Fig. [5] for sub-GeV neutrinos 
in angular bin —0.8 < cos?? < —0.6. The solid [black] 
curve employs our best fit parameters. To show the effect 
of 6'i3, we also plot the e-like events for 6>i3 = ±0.15 with 
A23 and 6*23 unchanged. The dash-dot [red] curve has 
negative 6*13, and the dashed [blue] curve has positive 
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1/E (GeV ) 



FIG. 6: [color online] The experimentally measured quantity 
TZ e versus the inverse neutrino energy E~ x in the sub-GeV 
region for angular bin —0.8 < cos •& < —0.6. The solid [black] 
curve is the result for the best fit parameters, the dash-dot 
[red] curve is for #13 = —0.15 and the dashed [blue] curve is 
the result for #13 = +0.15. 



6*13. It is clear that a negative value of this mixing angle 
permits an excess of e-like events for sub-GeV neutrinos. 

Returning to Fig. 01 we see how the full three-neutrino 
oscillation probabilities jointly affect the allowed region 
for the two mixing angles #23 and #13 . Using the subdom- 
inant approximation, the dash-dot-dot [blue] curve rep- 
resents the allowed region when all the data is employed. 
As expected, the region is symmetric about 6*13 = 0, and 
we note that the actual value of #13 has little impact upon 
the allowed values of #23, save the neighborhood immedi- 
ately around $13 = 0. Inclusion of the higher order terms 
in the oscillation probability dramatically alters this pic- 
ture. As discussed above, the data now favor negative 
9 13, with the atmospheric Super-K data shrinking the 
90% CL contour, the solid [black] curve, for positive #13. 
No longer is the contour symmetric about a particular 
value of 813; hence, the true value of this mixing angle 
will impact the allowed region for the O23 mixing angle. 
In particular, the allowed region for 623 shrinks as #13 ap- 
proaches positive values. In the future, should a reactor 
neutrino experiment confirm a nonzero value for |#i3|, it 
will have interesting consequences for the allowed value 
of 623- With such a measurement, we would perhaps 
see two true local minima in the A% 2 versus #13 plot in 
Fig. [31 The impact upon Fig. Q] would be to separate 
the jointly allowed regions into two disconnected curves 
with the limit on 623 more tightly constrained for positive 
values of 813. 



IV. CONCLUSIONS 

As we enter into the era of precision neutrino experi- 
ments, small effects, such as those arising from #13 or the 
octant of #23, require a careful treatment in the analysis. 



Future reactor experiments [29j, |3Cj, l3l| are sensitive to 
#13 and thus can determine the magnitude of #13, but not 
its sign. Long-baseline experiments, e.g., [3i| . will con- 
tain small effects that arc linear in #13, while an upgraded 
Super-K will produce additional data in the region where 
we have found significant effects linear in 613. How these 
different data interplay with each other in determining 
#13, including its sign, and the octant of 623 will be most 
interesting. 

In Ref. [49[ , it was shown that the constancy of R e im- 
poses an upper bound on 1 6^13 1 as well as constrains #23 
to be near maximal mixing. We find that present atmo- 
spheric data restrict the value of #13 from above, while 
the limit from below remains as determined by CHOOZ. 
We find 6>i3 = -0.07±^, assuming no CP violation. We 
have investigated which data points lead to the asymme- 
try in #13 and find that it is the atmospheric data in the 
very long baseline region previously noted [H], HfJ [HI to 
have significant terms linear in #13. We further found 
that the Earth's MSW effect plays an important role as 
it increases the effective value of #^2 m matter such that 
the atmospheric data provides a strict upper bound on 
#13. Further, the data producing the preference for a neg- 
ative #13 is data with an excess of e-like events, lZ e > 1. 
Allowing #13 to be negative supports this excess and per- 
mits #23 to be in the second octant. The parameters #13 
and #23 are found to be correlated; the statistically in- 
significant negative value for the minimum of #13 relates 
to the minimum for #23 being statistically insignificantly 
in the second octant, and the error in #23 is dependent on 
the value of #13 . Future measurements of #13 will impact 
the allowed value for #23. 

For A32 and #23 we find A32 
#23 - tt/4 = 0.03l£;^, where the use of the full three- 
neutrino oscillation probabilities leads to the asymmetry 
in the errors. We find that a quantitative analysis re- 
quires utilizing the more finely binned atmospheric data 
of Ref. [ill , the use of the full three-neutrino oscillation 
probabilities, and the inclusion of the full three-neutrino 
MSW effect. 



APPENDIX: EXPERIMENTAL SIMULATION 

In this appendix, we present the computational tools 
we use to analyze the Super-K atmospheric, CHOOZ, 
K2K, and MINOS experiments. The analysis tool for 
the Super-K atmospheric data is similar to that being 
used by others [HI, [2?], HI]; however, it is distinct in 
that we employ a full three-neutrino oscillation proba- 
bility rather than an approximate expansion, use a full 
three-neutrino treatment of the Earth's MSW effect, and 
include a model of the multi-ring data. The analysis of 
CHOOZ, K2K, and MINOS data is standard. Additional 
details can be found in Ref. L 50!]. Also, in this appendix, 
we demonstrate the efficacy of our analysis tools by com- 
paring our results with others' when appropriate. 

The appendix is organized as follows. We first discuss 
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the Super-K atmospheric experiment, beginning with the 
contained events followed by the upgpoing muon events. 
Then we discuss our statistical treatment of this exper- 
iment. Finally, we include a similar discussion for the 
CHOOZ, K2K, and MINOS experiments. 



1. Super-K contained events 

In order to observe atmospheric neutrinos at Super- 
Kamiokande, the neutrinos must interact with matter 
in either the detector or the surrounding environ to 
produced charged particles. The direction and energy 
of these charged particles can be deduced from the 
Cherenkov light they emit while traveling through the 
water- filled detector; from this data, one can infer on 
average the direction and energy of the initial neutrino. 
The Super-K experiment classifies the various detections 
in terms of the production point of the charged lepton, 
the number of charged particles produced, and their sub- 
sequent motion through the detector. 

Contained events refer to events in which the charged 
lepton is produced by the neutrino within the detector. 
These events are subdivided into fully contained and par- 
tially contained events. If an event is fully contained, 
then the charged lepton(s) produced within the detector 
do not escape the detector. An event is partially con- 
tained if the charged lepton(s) exit the detector. Finally, 
these two data sets are further separated into single-ring 
and multi-ring events according to the number of charged 
particles produced by the neutrino; if only one charged 



lepton is observed in the detector, it is termed a single- 
ring event. We first discuss the fully contained single-ring 
events and then extend this analysis to the other classes 
of data. The fully contained single-ring events are sta- 
tistically the most significant subset of the data and the 
cleanest to analyze. Preliminary discussions of the anal- 
ysis technique utilized for the fully contained events can 
be found in Refs. [5ll.l52j. 

The Super-K detector distinguishes between electrons 
and muons by the fuzziness of the Cherenkov ring gener- 
ated by the charged lepton; however, the detector cannot 
differentiate an electron e~ from a positron e + or a fi~ 
from a Given flavor and charge conservation, the de- 
tector can only determine if an event is e-like, originating 
from either a v e or v e interaction, or /i-like, originating 
from either a or F M interaction. Thus the detector 
counts charged leptons of flavor a in energy bin m and 
zenith angular bin n over the run time T 



nm 
a — >a 



at 



dt 



T 



(A.l) 



where the quantity dNp™ a /dt represents the rate at 
which a neutrino, created in the atmosphere with flavor 
/3, will be detected as an a-like event in the appropriate 
energy and angular bins within the detector. 

This rate depends upon the atmospheric neutrino flux, 
the neutrino oscillation probability from source to detec- 
tor, the kinematics of the charged lepton production, and 
the detector efficiencies. We may write it as 



— = N J dE ms J dcos9 u J dE v j dcos(9 s 

d 2 $ a (g„, cos 9 V ) - d 3 a p {E v ; Eg, cos6 s ) 

x £{E vis ) — — — P a /3(E U , cos9 v ) — — — — — 

dE v d cos 9 V dE e d cos 9 S d(j) s 

x e(cosiT< mQa; -cos??) e(cos7?-costf n - m<n ). (A.2) 

I 



We define the variables in Eq. (|A.2|) . N represents the 
number of target protons. E V i S is the energy measured 
by the detector (this quantity is defined differently de- 
pending on the data sample); E™'J nax (E™' s mm ) is the 
maximum (minimum) value of E V i s for bin m. (For 
single-ring fully contained events, E V i S is simply the en- 
ergy of the created lepton, E V i s = Eg.) is the zenith 
angle of the detected charged lepton with cos i3 = 1 in- 
dicating the vertically downward direction. The rela- 
tive angle between the incident neutrino and the pro- 
duced charged lepton are described by the 9 S scattering 
angle and the cj> s azimuthal angle. The energy of the 
incident neutrino is E v with zenith angle 9 V . The az- 
imuthally averaged atmospheric neutrino flux for a neu- 



trino of flavor a is d 2 $ a (E v , cos9 v )/dE u dcos 9 V which 
we take from Ref. [53]. e(E V i S ) corresponds to the de- 
tection efficiency. & 3 ap(E v ; Eg, cos0 s )/dE^ dcos# s d<^ 
is the differential cross section for a neutrino of energy 
E v and flavor (3 to produce a charged lepton of flavor j3 
with energy Eg through a scattering angle 9 8 . Although 
the differential cross section which occurs in Eq. lA.2l does 
not depend on the azimuthal angle <p s , the geometry that 
determines in which angular bin an event lies does depend 
on <j> a . This is because the zenith angle ■& of the charged 
lepton is given in terms of the neutrino zenith angle 9 V 
and the scattering angles 9 S and <f> a by 

cos i? = cos# s cos 9 U — sm9 s sin6*„ cos</> s . (A. 3) 
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The energy range for atmospheric neutrinos as mea- 
sured at Super-K requires the use of several cross sec- 
tions. At low energies, below 1 GeV, the dominant 
process is charged-current quasi-elastic scattering from 
the proton and the nucleons in the oxygen nucleus in 
H2O, e.g., V e + p + — > e + + n and v e + n — > e~ + p + . 
At intermediate energies, peaking around 1.5 GeV, the 
dominant process is single-pion resonance production, i.e. 
v a +N -> £ a +N* followed by N* ->■ N'+tt. At higher en- 
ergies, starting at 1 GeV and dominating above 10 GeV, 
deeply inelastic scattering occurs, v a + N — > £ a + X 
We utilize the same set of cross sections as was used 
111 Ref. 0. The Heavisidc functions are inserted into 
Eq. (|A.2p in order to restrict the values of 9 V , S , and (j) s 
to values which produce a value for 1? that lies within bin 
n. 

The incident neutrino's zenith angle 9 V does not 
uniquely determine the path length L, as neutrinos are 
produced at a variety of vertical heights h in the at- 
mosphere. We thus introduce the oscillation probability 
Vap{E v , cos 9 V ) averaged over this production height in 
the atmosphere 

P a /3(E V , cos 9 V ) = 

/ dhP a (h,E„)V a0 (L(h, cos 9 V )/E V ), (A.4) 
Jo 

where L is related to h and cos 9 V by 

L = y/R 2 cos 2 9 V + h (2R + h) - R cos9 u , (A.5) 

with R the radius of the Earth. P a (h, E v ) is the normal- 
ized probability for a neutrino of flavor a to be created 
at a height h, a quantity we take from Ref. [H3 |. 

The neutrino oscillation probability in vacuum 
~Pap(L/E v ) is given in Eq. (JTJ; however, the coherent 
forward scattering of neutrinos on matter alters the prob- 
ability for those neutrinos which pass through the Earth 
[4TI [43 ] . Neutral current interactions between the neutri- 
nos and matter are not flavor dependent leaving the oscil- 
lation probabilities unaffected; however, charged current 
interactions will introduce into the Hamiltonian a flavor 
dependent potential. In the flavor basis, we may write 
the neutrino evolution equation as 



where M = diag (0, A 2 i, A31) and the operator V op- 
erates on only the electron flavor with a magnitude 
V = \/2GfN,,(x) in which N e (x) is the electron number 
density of the matter. We note that for anti-neutrinos, 
this potential has the opposite algebraic sign. In matter 
of constant density, one may diagonalize the Hamilto- 
nian to determine the effective mass-squared differences 
in a new propagation basis; this basis is related to the 
flavor states by a modified effective mixing matrix in 
the matter. The density of the earth may be approxi- 
mated as piecewise constant [43[ so that neutrinos which 
travel through the earth may traverse regions of varying 
density. For certain energies and mixing angles, neu- 
trinos can undergo parametric resonances when passing 
through regions of different densities [44|. For a piece- 
wise constant density profile, we may use the the meth- 
ods of Refs. [H, H3| to exactly determine the neutrino 
oscillation probability. This method is computationally 
efficient as it merely involves finding the local effective 
mass-squared differences and mixing angles, and its ex- 
actness accurately accounts for effects which depend on 
small parameters or effects due to parametric resonances. 



id t Vf 



1 

2E~u 



UMU^ + V 



(A.6) 



The full calculation represented by Eq. (|A.2[) is a nu- 
merically intensive five-dimensional integration. We de- 
sire an analysis tool which is both numerically accurate 
and sufficiently computationally efficient that we can 
scan a large swath of parameter space in a reasonable 
amount of time. To effect this, we make approximations 
in Eq. (|A.2[) regarding the scattering angle and efficiency 
terms. First, we fix the scattering angle 9 S to its average 
value 9 S as a function of E V i S = Eg for each energy and 
angular bin as calculated by the Monte Carlo simulations 
in Ref. [56j . This eliminates the integral over cos 9 S . We 
are justified in doing so because, if this integral is done 
last, the integrand is quite smooth and nearly linear over 
each bin, particularly for the finer binning of the most 
recent data. Thus, we expect this approximation to be 
quantitatively accurate, and this is what we have found. 
With this approximation, Eq. (|A.2[) becomes 



f^]\[ nm poo pE ' p 

= N J AEv j mmin dEmS J dC ° S6 

.„ , d 2 $ a {E v> cos 9 V ) ~ d 2 ap{E v \ E t ,cos6 3 (E via )) 

x £{E ms ) — — — V a p[E v ,cosO v ) — 

AE V d cos 9 V dE e d(j> s 

x e(cos?T< mQ;i; -costf) e(costf-costT' mm ), (A.7) 

I 



where the differential cross section and the angle d have been evaluated at cos(# s ) = cos(0 s ). 
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Finally, we must determine the detector efficiency provided. For no oscillations, Eq. (|A.2[) becomes, utiliz- 
e{E vis ) 1 which has not been furnished by the experimcn- ing our assumptions, 
talists. We can, however, extract it from information 



dN™ f°° [ E ?u max r f dj> a (^,cos^) d 3 a [3(E„; E e , cos6 s (E ms )) 

- e[m)N \ dE v I dE ms / dcos6>„ / d<j> s — — ^— — ' , — -j — (A. 



dt J J E m,min J J dE u d cos U dEf d cos 



where we have assumed that the efficiency is primarily 
dependent on the lepton energy, e(E V i S ). Further taking 
it as a constant e(m) over each energy bin to, we pull 
it out of the integral. The Monte Carlo calculation of 
dN™/dt appears in Ref. [57|- Performing the integrals 
on the right hand side of Eq. (|A.8[) , we can use the Monte 
Carlo result to determine the average efficiency e(m) for 
each energy bin. This average is then used in Eq. (|A.7|) . 

To conclude the discussion of the single-ring fully con- 
tained analysis, we note that there are ten angular bins 
of equal size in cos?9, going from +1 (downward) to — 1 
(upward); ten energy bins for e-like events; and 8 energy 
bins for /z-like events. In total, this subset of the data 
consists of 180 data points. 

Within the fully contained data set, we also have multi- 
ring events. This data correspond to neutrinos which 
interact inside the detector to produce more than one 
detected particle in the final state. We can calculate the 
multi-ring event rate in a manner similar to the single- 
ring event rate provided we modify the visible energy and 
neutrino event direction. In a multi-ring event, the visi- 
ble energy depends on the number of particles in the final 
state, their momenta, and their scattering angles. There 
is not a simple and reliable way of determining E V i S ; as 
such, we must make some approximations. We use the 
results of Monte Carlo simulations in Ref. (5(| to esti- 
mate the average value of E v / E V i S for each energy bin. 
To determine the angular distribution of these events, we 
assign a single scattering angle to the final state particles. 
To make a simple estimate of the average multi-ring scat- 
tering angle for each energy bin m, 9 Stmu iti(rn), we fit our 
no-oscillation event rate to the analogous no-oscillation 
Monte Carlo calculation from Ref. [57| . Our analysis then 
uses Eqs. (|A.llA.2IA.3p to calculate the general event 
rate. 

The next data set consists of the partially-contained 
events. In these, the charged lepton is created within 
the detector but has sufficient energy to escape the inner 
detector and be detected by the outer detector (OD). If 
the energy deposited in the OD appears to be less than 



that needed for a muon to traverse the OD, the event 
is classified as an OD stopped event; those depositing 
more energy are termed OD through-going. We alter the 
definition of the visible energy 

Eyis Ei nner -\- Edead E an n , (A. 9) 

with Ei nner the total energy of the charged particle ob- 
served in the inner detector, E^ ea d the energy deposited 
in the region between the inner detector and the outer 
detector, and E ant i the energy observed in the outer de- 
tector. The charged lepton energy is then a function of 
this visible energy; we set Eg{E V i S ) to its average value as 
taken from the Monte Carlo calculations from Ref. [5t| . 
The efficiencies are similarly taken for each data bin from 
the same Monte Carlo no-oscillation calculations as de- 
scribed above. For these single-ring partially contained 
events, the same ten angular bins are used, and there 
are four energy bins each for both the OD stopped and 
through-going events, giving a total of 80 data points. 



2. Upgoing muons 

A second major subset of the atmospheric data arises 
from detecting muons created by atmospheric neutrinos 
and anti-neutrinos interacting with the rock surrounding 
the detector. These muons can then either stop in the 
detector, called "stopping muons" , or pass through the 
detector, called "through-going muons" . Typically the 
stopping muons have an O(10 GeV) energy while the 
through-going muons have an 0(100 GeV) energy. Al- 
though statistically not as significant as the contained 
events, these data provide important information in this 
higher energy range. 

Our analysis of these events exactly follows that given 
in Refs. [58L Hl|. We do not include the effects of the 
muon energy fluctuations as has recently been done in 
Ref. [27|. The detection rate for stopping muons, S, and 
through-going muons, T, is given by 
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dN 



H,S,T 

dt 



d cos 6,, 



d 2 $0,(^,008^) - 



dE v d cos t 



V a ^(E v , cos6>„) 



da^(E v ; 
dE„ 



TZ(E^E th ) A s , T {E t h, cosOy) . 



(A.10) 



We have assumed, appropriate for these high energies, 
that the scattering is forward. This allows us to replace 
the charged lepton angle d with the neutrino angle B v \ 
we may also perform the integration over the scatter- 
ing angles O s and <j> s in the cross section. Na is Avo- 
gadro's number. The function 7Z(E^,,Eth) is the aver- 
age distance that a muon of energy E^ will travel un- 
til its energy reaches the value E t h, the amount of en- 
ergy needed to traverse the detector; this quantity is ex- 
pressed in the natural units for range, distance times 
the Earth's density. At (Eth, cos 8 V ) is the area pro- 
jected onto a plane perpendicular to the muon direc- 
tion such that a muon of energy E t h or greater can pass 
through this part of the detector. The details for cal- 
culating AT(E t h,cosO v ) can be found in Ref. [60| . For 
the stopping muons, As(E th ,cosO u ) = A(£ min ,cosd v ) - 
At (Eth, cos 8 v), where A(£ m i n ,cos& v ) is the projected 
area of the detector with a path length greater than i m i n 
taken to be 7 m by the experimentalists. Note that there 
is only muon data as electrons/positrons produced in the 
rock are unable to travel to the detector. The data cov- 
ers the angular region from cos#„ = to cos 9 U = —1, 
directions where muon production from the rock exceeds 
the cosmic ray background. Since the neutrinos can orig- 
inate as either electron or muon neutrinos, we sum over 
the two neutrino flavors a as in Eq. (|A.10|) ■ The upgoing 
muon data is binned in ten angular bins and not binned 
in energy, resulting in a total of 20 data points. 

Analysis of these events is not as computationally in- 
tensive as the calculation of the contained events because 
the forward scattering allows the integration over the 
scattering angles and the muon energy E^ to be per- 
formed outside the fitting program. The parameters be- 
ing fit are contained in V atl (E v ,cos6 y ) which is inde- 
pendent of the muon energy and scattering angle. We 
have found it efficient to change the integration over these 
variables to the Feynman scaling variables x and y, as is 
natural for the deeply inelastic region. 



3. Super-K Statistical Analysis 

We use the most recent experimental data from 
Ref. [l6| which includes 180 data points for fully con- 
tained single-ring events, 90 for fully contained multi- 
ring events, 80 for partially-contained events, 10 for up- 
ward through-going muons, and 10 for upward stopping 
muons; in all, this constitutes 370 data points. In order 
to determine the neutrino oscillation parameters, we con- 




FIG. 7: Allowed region at 90% confidence level for A32 and 
sin 2 2023 in the subdominant approximation using only Super- 
K atmospheric data. 



struct a x 2 based upon a Poisson distribution, following 
the same procedure used in Ref. [l6j]. We incorporate 
systematic errors by utilizing the "pull" approach as de- 
scribed in Ref. [? ], which allows one to incorporate sys- 
tematic errors in the analysis without adding adjustable 
parameters. The approach is based upon allowing linear 
corrections to the theoretical predictions for each system- 
atic error. Our \ 2 function is 



:i7<> 



E 

n=l 

+2N obs (n) In 



2{N the (n) - N obs (n)} 

N obs (n) 
N th e(n) 



43 ,m2 

l — l X 



(ATI) 



Nobs(n) is the number of observed events in the bin n; 
Nthe( n ) is the theoretical prediction of the number of 
events in that bin; £i is the systematic error pull for the 
systematic error i; and ui is the one-sigma value for the 
systematic error i. Nthe(n) represents a modified pre- 
diction of the expected number of events due to the in- 
clusion of systematic errors; the systematic errors adjust 
this quantity through an assumed linear dependence on 
the pulls £i. Here we use 45 systematic errors arising 
from different inputs into the data analysis as described 
in Tables VII-X taken from Ref. (l5| . For these 45 er- 
rors, all of them contributed to the \ 2 except the overall 
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flux normalization and the normalization for the multi- 
GeV multi-ring sample, which are floated freely. During 
each fit, these forty- five & are varied to minimize x 2 for a 
given set of oscillation parameters. The minimization of 

X 2 with respect to & (^- = 0) is equivalent in the pull 
method to numerically solving for & in the 45 coupled 
equations of the form 



400 - 



dx 2 




(n)- 



0, 



N bs (n) 



(A.12) 



where /" is the zth systematic error for bin n. These 
equations are linearized and solved as a set of coupled 
linear equations. 
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FIG. 8: [color online] Ax 2 versus sin 2 2#23 for our analysis 
of the recent Super-K atmospheric data in the sub- dominant 
approximation. The horizontal lines take the values: [red] 
solid, Ax 2 = 2.30; [blue] dashed, Ax 2 = 4.61; and [green] 
dot-dashed, Ax 2 = 9.21, the 68%, 90%, and 99% confidence 
levels for a two parameter fit. Both A32 and #13 are varied in 
obtaining this curve. 
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FIG. 9: The top two graphs depict the number of events for 
the sub-GeV data as a function of the angular bin in cosi?, 
the middle two graphs are the same except for the multi-GeV 
data, and the bottom two graphs are for the multi-ring events. 
The [black] solid curves represent the results from our best fit 
parameters, the [red] dashed curves represent the data, and 
the [blue] dot-dashed curves represent the Monte Carlo no- 
oscillation predictions. 



solid curve in Fig. 3a, —0.38 < 6*13 < 0.38. This is exactly 
the result of Ref. [lj| , and our other results are in excel- 
lent agreement with that analysis. As noted previously, 
our reproduction of the allowed region for #13, which has 
a nonzero but small effect on the atmospheric data, is a 
very strong test of our analysis. 

Finally, in Figs. O and [TU1 we compare the predicted 
number of neutrino events corresponding to our best fit 
parameters with the experimental data as a function of 
the zenith angle. We also present the Monte Carlo pre- 
dictions for the expected number of events in the absence 
of neutrino oscillations. Each of the different Super-K at- 
mospheric data sets is depicted. The results are a good 
fit to the data, comparable to that found in Ref. jlq . 



We compare our analysis with that performed by the 
experimentalists in Ref. fl6j |. To do this, we utilize the 
sub-dominant approximation in Eq. ([2]) and minimize 
the above x 2 - Our best fit oscillation parameters are 
(A 32 , sin 2 6*23,6*13) = (2.5 x 10- 3 eV 2 , 0.51, 0.01) with an 
overall x 2 of 416 for the 370 data points. In Fig. l.a, the 
[black] solid curve represents Ax 2 versus A32 for the sub- 
dominant approximation, using only atmospheric data. 
In Fig. [Jj we present the allowed region for A32 and 
sin 2 26*23 corresponding to Ax 2 = 4.61, the 90% confi- 
dence level for a two parameter fit. We also present Ax 2 
versus sin 2 6*23 in Fig. [5) 

At the 90% confidence level for a two parameter fit, 
we find the allowed parameter values 2.1 x 10~ 3 eV 2 < 
A 32 < 3.1 x 10" 3 eV 2 and 0.938 < sin 2 26» 23 . Addition- 
ally, we extract the allowed value for 6^3 from the [black] 



4. CHOOZ experiment 

For the CHOOZ reactor experiment, we follow a stan- 
dard procedure as described in Ref. [211 ]. In our analysis, 
we use experimental data that consists of seven positron 
energy bins for each of the two reactors, giving a to- 
tal of 14 bins. We include a 14 x 14 covariance matrix, 
VjJ , to account for the correlation between the energy 
bins, and we include the systematic error from the over- 
all normalization and energy calibration. We write the 
expected positron yield for the k th reactor and the j th 
energy spectrum bin as 

X{E 3 ,L k ,e, A 32 ) = XiEjjPiE^L*, e, A 32 ), 

(j = l,...,7, k = l,2), (A.13) 
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where X(Ej) is the distance-independent positron yield 
in the absence of neutrino oscillations, is the reactor- 
detector distance, and P{Ej, Lk, 0, A32) is the oscillation 
probability averaged over the energy bin and the detector 



and reactor core sizes. In our fitting routine, we minimize 
the following \ 2 function with respect to the neutrino 
oscillation parameters 



X 2 (0,A 32 ,a,g) = 

14 14 

J2 ( X i " a ~X(9Eu U, 9, A 32 )) V (Xj - aX(gEj, L 3 , 6, A 32 ] 

i=l j=l 



1 \ 2 / 1 \ 2 

a — 1 \ 1 g — 1 



, (A.14) 



with the absolute normalization constant a and the 
energy-scale calibration factor g. 
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FIG. 10: The same as Fig. [9]except the data sets are now the 
partially contained stopping events (upper left), the partially 
contained through-going events (upper right), the upward 
stopping muon events (lower left), and the upward through- 
going muon events (lower right). 



5. K2K experiment 

For the K2K experiment [20], we employ the method 
developed in Ref. [6l| to estimate the expected no- 
oscillation neutrino spectrum, S{E V ) 1 in the relevant en- 
ergy range of ~ 0.2 to ~ 3.0 GeV. The expected number 
of neutrino events for oscillating neutrinos is then 



Xf eo = / SiEJP^L/Ev) 

JE min (n) 



(A.15) 



where P^^L/Ev) is the muon neutrino survival probabil- 
ity and E max (n) (E m i n (n)) are the maximum (minimum) 
energy values for the energy bin n. For the statistical 
analysis, we follow the procedure described in Ref. 
We only use the single-ring sub-sample, which consist of 




, (10" 3 eV 2 ) 



FIG. 11: [color online] A\ 2 versus A32 for our analysis of the 
K2K experiment in the sub-dominant approximation. The 
horizontal lines take the values: [red] solid, Ax 2 = 2.30; [blue] 
dashed, A\ 2 = 4.61; and [green] dot-dashed, A\ 2 = 9.21, the 
68%, 90%, and 99% confidence levels for a two parameter fit. 



58 neutrino events. The signature for neutrino oscilla- 
tions from i/fi to v T in a two neutrino analysis are both a 
reduction in the total number of observed neutrino events 
and a distortion in the neutrino energy spectrum. The 
X 2 function is divided into two terms: the observed to- 
tal number of events detected at the Super-K detector, 
Xnormi an d the shape of the spectrum included in xlhape- 
We use the "pull" method [? ] to account for 31 system- 



atic uncertainties by adding a third term x 



syst ' 



2 _ 2 
XK2K — Xnorm 



Xshapc 



Xsyst 



(A.16) 



The best fit oscillation parameters, A32 and 623, are ob- 
tained by minimizing Xk2K- 

The systematic parameters included in x 2 ys t arise from 
the neutrino energy spectrum at the near detector site, 
the flux ratio, the neutrino-nucleus cross-section, the ef- 
ficiency and the energy scale of the Super-K detector, 
and the overall normalization. The kth. systematic error 
is represented by the coefficient C% and modifies the ex- 
pected number of neutrino events, Eq. (IA.15[) . in a linear 
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manner according to the "pull" method 

31 



3.5 



N, 



thco 



thco 
total 



fc=i 



with £fc the pull corresponding to systematic error k. 

Due to the low statistics, we employ a Poisson distri- 
bution; hence, the expressions for Xnorm an( i X 
given by 



shape 



are 



Xnorm 



Xshapc 



2 N, 



thco 
total 



data 
total 



n: 



data 
total 



/ythco 

ln% 

/vrdata 
Jv total 




data 



iVf ata ln 



n; 



the 



(A.18) 



where A^ ata is the exp erimental data provided by the 
K2K collaboration [20J and the superscript "total" im- 
plies a sum over n. The contribution to y 2 from the 
systematic errors is 



Xsyst 



31 
j,k=l 



(A.19) 



where we use an error matrix Mk i co nstructed from Ta- 
bles 8.1 and 8.2 provided in Ref. 

In Fig. [TTJ we depict A% 2 versus A 32 for an analy- 
sis that utilizes only the K2K data in the subdominant 
approximation; also, Fig. [T2] shows the 90% CL allowed 
region in the [red] dashed contour for A32 and sin 2 2# 23 in 
the subdominant approximation. The absolute minimum 
in our fit is (A 32 , sin 2 26> 23 ) = (2.78 xl(T 3 eV 2 , 0.998). At 
90% CL, we find 2.2 x 10~ 3 cV 2 < A 32 < 3.2 x 10~ 3 eV 2 . 
The total number of observed events 58 is in agreement 
with the 56 events found from the model. All of these 
results are consistent with the analysis performed by the 
experimentalists in Ref. [20] |. 



6. MINOS experiment 

The MINOS experiment is quite similar to the K2K ex- 
periment and thus we apply a similar analysis technique. 
For the no-oscillation spectrum, we use the Monte Carlo 
simulation provided by the MINOS collaboration [lg| . 
We normalized this spectrum to measurements made at 
the near detector. A total of 1065 events were expected 
in the absence of neutrino oscillations. With the no- 
oscillation spectrum and Eq. (|A.15[) , we calculate the 
expected number of neutrino events in the presence of 
neutrino oscillations. We use the MINOS data [l8| cor- 
responding to two years of beam operation in which 884 
neutrino events are observed. The data consists of 
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FIG. 12: [color online] A32 vs sin 2 2#23 allowed region for a 
two parameter fit in the subdominant approximation at the 
90% CL for the K2K data, [red] dashed, and the MINOS, 
[blue] solid, experiments. 




FIG. 13: [color online] The same as Fig. 1111 except the data 
are from MINOS. 



15 energy bins along with three systematic errors: the 
relative normalization between the far and near detec- 
tors with a 4% uncertainty; the absolute hadronic energy 
scale with a 11% uncertainty; and a 50% uncertainty in 
the neutral-current background rate. For the definition 
of x 2 ) we use the Poisson distribution function 
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thco ^y-data ^y-data n 




Xminos = 2 E ( ^ 

n=l 
3 

+E 

i=i 



where the symbols are analogously defined to those in 
the K2K section. 

In Fig. rj3] we plot Ax 2 versus A 32 for the K2K data 
in the subdominant approximation; likewise, in Fig. 1121 
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we present the allowed region at 90% CL, the [blue] solid 
curve, for A 32 and sin 2 2# 2 3- The minimum of x 2 is lo- 
cated at A 32 = 2.41 x 1(T 3 eV 2 and sin 2 26» 23 = 0.99 90. 
The allowed intervals of these parameters at 90 % CL 
are 2.25 x 10~ 3 eV 2 < A 32 < 2.8 x 10~ 3 eV 2 and 
0.86 < sin 2 26>23 for A^ 2 = 4.6. All of these results are 
consistent with the analysis performed by the experimen- 
talists in Ref. [H|. 

Results which combine CHOOZ, K2K, and MINOS are 
presented in the main body of this work. This Appendix 
provides details of the analysis tools we use throughout 
this work; additional details may be found in Ref. [501 ] . 
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